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Equation of state of the unitary Fermi gas: An update on lattice calculations 
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The thermodynamic properties of the unitary Fermi gas (UFG) have recently been measured 
to unprecedented accuracy at the MIT. In particular, these measurements provide an improved 
understanding of the regime below T/cp ~ 0.20, where a transition into a superfluid phase occurs. 
In light of this development, we present an overview of state-of-the-art auxiliary field quantum Monte 
Carlo (AFQMC) results for the UFG at finite temperature and compare them with the MIT data 
for the energy, chemical potential, and density. These AFQMC results have been obtained using 
methods based on the hybrid Monte Carlo (HMC) algorithm, which was first introduced within the 
context of lattice QCD. 
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The unitary Fermi gas (UFG) is defined as a two- 
component many-fermion system in the limit of short 
interaction range Tq and large s-wave scattering length a, 
such that ^ fep^^n <^ 1 <SC k^a — ^ 00, with 



kp = (Stt^tt,)^/^ being the Fermi momentum and n be- 



ing the particle number density (we choose units such 
that h = kg — 1). The UFG also saturates the unitarity 
bound on the quantum-mechanical scattering cross sec- 
tion (To(fc) < 47r/fc^, where k is the relative momentum of 
the colliding particles. The UFG features special prop- 
erties that arise from the fact that it is characterized by 
a single length scale, given by the interparticle distance 
^ kp^, independently of the details of the interaction. 
While the thermodynamic properties of the UFG are 
universal the lack of a readily accessible dimension- 
less expansion parameter renders the UFG a challenging 
many-body problem. Since the proposal of the UFG as 
a model for dilute neutron matter by Bertsch [2| and its 
realization in ultra-cold atom experiments (see Ref. Q 
for a review of the experimental situation), the UFG has 
received widespread attention across multiple disciplines, 
ranging from atomic physics p| to the study of nuclear 
matter Q and relativistic heavy- ion collisions [6] . 

On the experimental side, the presence of a super- 
fluid phase in the UFG below T/ep ~ 0.15 was demon- 
strated directly a few years ago through the creation of an 
Abrikosov vortex lattice under rotation [7]. However, a 
direct thermodynamic signature of the transition was not 
unambiguously established until the recent high-precision 
measurement at MIT of the equation of state (EoS) of 
a homogeneous two-component UFG over a large tem- 
perature range Q . These measurements were performed 



on trapped ^Li atoms (using a Feshbach resonance to 
tune the system to the unitary limit), which enabled a 
detailed study of the compressibility, density, and pres- 
sure of the UFG. In addition, greatly refined empirical 
results were obtained for the associated critical temper- 
ature T^/ep = 0.167(15) as well as for the "Bertsch pa- 
rameter" ^ ~ 0.376(5), which characterizes the ground 
state of the UFG. As precision data are now available for 
the energy, chemical potential, and density of the UFG 
in a wide temperature range, an opportunity presents it- 
self to compare these measurements with calculations in 
various theoretical frameworks. Here, we focus on com- 
paring the MIT data with the most recent auxiliary field 
quantum Monte Carlo (AFQMC) results. 

The Hamiltonian that captures the physics of the uni- 
tary limit can be written on a spatial lattice as 



(1) 



where A denotes the spin projection, m is the fermion 
mass (we also set m = 1), and g is the coupling constant. 
The creation and annihilation operators satisfy fermionic 
anticommutation relations, and fixi'^'i) = '^xi'^'i) 
denotes the number density operator at lattice position 
rj for spin projection A. The thermodynamic equilibrium 
properties are obtained from the partition function 



Z = Tr exp[-(3{H-fiN)] 



(2) 
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where N is the total particle number operator, n is the 
chemical potential, and /? = 1/T is the inverse tempera- 
ture. 

To evaluate expectation values of observables numeri- 
cally, we followed the path-integral approach presented 
extensively in Ref. 0], with recent improvements de- 
scribed in Ref. [l3|- The system is placed on a cubic 
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spatial lattice of extent L = N^l with periodic bound- 
ary conditions. The lattice spacing / (henceforth set to 
unity) and extent L provide natural ultraviolet (UV) and 
infrared (IR) momentum cutoffs, given by k^^^ = it /I 
and /cq — 27r/L, respectively. The imaginary-time evo- 
lution operator exp[— /3(if — ^N)] is expanded using a 
Trotter-Suzuki decomposition with temporal lattice spac- 
ing T, and the interaction is represented by means of a 
Hubbard- Stratonovich (HS) transformation [ll|. As we 



focus on the spin-symmetric case, the fermion sign prob- 
lem is absent. The resulting path integral formulation is 
an exact representation of Eq. ([2|) up to finite- volume and 
discretization effects, which may be controlled by vary- 
ing the spatial lattice volume V = N!^ and density n. 
The thermodynamic and continuum limits are recovered 
as y — >■ cxD and n — >■ 0, respectively. The latter requires 
great care, as too low densities may introduce shell ef- 
fects. As our lattice formulation is very similar to that of 
Ref. referred to here as determinantal Monte Carlo 
(DMC), we shall restrict ourselves to describing three 
modifications which significantly improve the results. 

First, the bare-lattice coupling constant g correspond- 
ing to the unitary regime is determined by means of 
Liischer's formula [l^ as in Ref. [l^. This procedure 
yields g ~ 5.14 in the unitary limit. Our lattice Hamil- 
tonian contains g as the sole parameter characterizing 
the interaction. Finite-range effects are induced by the 
presence of the UV cutoff of the lattice. In order to min- 
imize such discretization effects, the dilute limit should 
be approached as closely as possible. Recent theoretical 
developments 14 - 1^ have explored the use of improved 
transfer matrices and operators, with multiple param- 
eters tuned to unitarity. The implementation of such 
methods is an objective of future AFQMC calculations. 

Second, we use a compact, continuous HS transforma- 
tion referred to as "Type 4" in Ref. 17|, which found 
it superior with respect to acceptance rate, decorrelation 
and signal-to-noise properties than the more conventional 
unbounded and discrete HS transformations fl8| . 

Third, we update the HS auxiliary field a using hybrid 
Monte Carlo (HMC), which combines molecular dynam- 
ics (MD) evolution of cr with a Metropolis accept-reject 
step [l^. The result is the determinantal hybrid Monte 
Carlo (DHMC) algorithm, introduced in Ref. jlO]. In 
DHMC, global MD updates of a take place via introduc- 
tion of a momentum field tt conjugate to cr, such that the 
dynamics is given by the Hamiltonian 
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(3) 



where h{[a] encodes the dynamics of the fermion degrees 
of freedom (for more details on U, see e.g. Ref. 0). The 
DHMC algorithm produces greatly enhanced decorrela- 
tion between successive MC samples for all temperatures 
and lattice sizes, and removes the necessity to spend an 
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FIG. 1: (Color online) AFQMC results for (top) the total 
energy E in units of the energy of a free Fermi gas and (bot- 
tom) the chemical potential fi in units of the Fermi energy 
as a function of T/ep. Results are shown for A^^. = 8 (green 
squares), A'^^, = 10 (red circles), = 12 (purple triangles), 
= 14 (blue inverted triangles), and = 16 (black aster- 
isks), where V = is the (spatial) lattice volume. 



increasing number of decorrelation steps at larger V (note 
that the computational cost of a full sweep of the lattice 
scales as ~ y in a local algorithm such as DMC). This 
is replaced by a fixed number of operations, typically of 
0(10^), required to produce one MD "trajectory", inde- 
pendently of V . 

The DHMC algorithm allows for an extension of the 
AFQMC analysis beyond the capabilities of DMC, which 
is currently limited to ~ 10 in the spatial lattice 
extent, n ~ 0.1 in particle number density, and N ~ 100 
in the particle number. The improved scaling of the CPU 
time has allowed us to study lattices up to = 16, while 
simultaneously maintaining a relatively large number of 
particles, N ~ 45,75,110, and 160 for = 10,12,14 
and 16, respectively, which corresponds to densities in 
the range n ~ 0.040- 0.045. The A^^ = 8 data (which are 
not identical to Ref. ^9]), correspond to ~ 35 and n ~ 
0.070, which was not reduced further in order to avoid 
shell effects. We generated ~ 200 uncorrelated snapshots 
of (7 for each value of T/ep, which yields a statistical 
uncertainly of ci; 1% for the observables. Expressions for 
the AFQMC computation of the latter were obtained by 
differentiating Z with respect to j3 and /i as conventional 
in thermodynamics, with the exception that Z is replaced 
by its discretized form in terms of the HS field. 

In Fig. [U we present AFQMC results (for various lat- 
tice sizes) for the total energy E in units of the energy 
of a free Fermi gas i?pQ = 3/5 A^e^, and for the chem- 
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FIG. 2: (Color online) Energy E/EpQ (red dots), as obtained 
by Ku et al. • Our AFQMC results extrapolated to infinite 
volume are shown by open black circles. The results for N,j. — 
8 (open blue squares) were obtained with the DMC algorithm 
in Ref. ^ . The green square shows the QMC result of Ref [2^ 
for ^ at T = 0. The inset shows the vicinity of the superfluid 
phase transition at T^/ep ~ 0.15. 



ical potential /i in units of €p = kp/2m as a function 
of T/ep. Before comparison with the MIT data we 
performed an extrapolation to the infinite volume limit 
Nx —7- oo. This required interpolation of the data series 
for each value of N,^^ as the physical temperature T/ep 
is not known beforehand. Apart from this minor com- 
plication, the extrapolation to infinite volume is greatly 
facilitated by the lack of a systematic variation in the 
results with the lattice volume above ~ 10. Since we 
have not performed an extrapolation to the continuum 
limit (which requires n — 0), our results may be affected 
to some degree by systematic errors due to the effective 
range r^f^. Our results currently reach kpr^^ ~ 0.3 which 
is non-negligible and may produce significant deviations, 
in particular at low T/ep as shown by Carlson et al. [2^. 

In Fig. [21 we compare our AFQMC results (extrap- 
olated to infinite volume) with the measured energy E 
of the homogeneous UFG. The overall agreement is sat- 
isfactory throughout the range of temperatures studied. 
At low T/ep, AFQMC slightly overpredicts the experi- 
mental data. Our new results show a noticeable improve- 
ment over the results of Ref. Q with = 8, likely due 
to decreased finite- volume and effective-range effects. In 
contrast to the case of E/E-pQ, our results for fi/ep in 
Fig. [3] deviate noticeably from experiment at low T/ep, 
where ^i/ep ~ 0.38 at T/ep ~ 0.1. AFQMC overpre- 
dicts this by ~ 5%, which clearly exceeds the statisti- 
cal uncertainty. However, the larger lattices used here 
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FIG. 3: (Color online) Chemical potential fi in units of ep 
as measured by Ku et al. [1]. The notation for the AFQMC 
results is identical to Fig. [2] The solid green square shows the 
result of Ref. 20'] assuming jj./ep{T = 0) = ^. 



represent a dramatic improvement over Ref. in par- 
ticular above T^/ep ~ 0.15. Nevertheless, the discrep- 
ancy below T^/ep cannot be accounted for at present. In 
Fig. m we show the particle number density relative to 
the temperature-dependent density of the noninteracting 
Fermi gas. Again, a discrepancy at low T/ep is found, 
which is analogous to that observed for [i/ ep. 

While the agreement between our AFQMC calculation 
and the data of Ref. Q is satisfactory in general, notable 
discrepancies persist. We have achieved a significant re- 
duction of the density from n ~ 0.1 to n ~ 0.04, with 
a concomitant decrease in discretization (finite-range) 
effects. Nevertheless, since finite-range effects scale as 
~ n^l'^ , this still only implies an effective reduction from 
n^/^ ~ 0.46 to ~ 0.34. The possibility that the discrep- 
ancies between our AFQMC data and experiment are due 
to residual finite-range effects can therefore not be ruled 
out at present. 

As the region where the discrepancies are largest ap- 
pears to be at very low T/ep (at least for E/E^q and 
/i/e^.), the task of performing calculations at significantly 
lower values of v}/'^ for such temperatures is extremely 
demanding, indeed largely beyond the capabilities of ex- 
tant algorithms. In this situation, accounting for the 
finite-range effects by improving the transfer matrix (as 



in Refs. [1J,|15|) provides a systematic way to remove the 



finite-range effects from both the action and the observ- 
ables at a given density, without modifying the temper- 
ature scale of the calculation. Preliminary results have 
appeared in Ref. [l6j . Another source of error under in- 
vestigation is the Trotter-Suzuki step r. This was found 
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FIG. 4: (Color online) Density n{ti, T) of the UFG (red cir- 
cles) as obtained by Ku et al. normalized to the density 
no{fj.,T) of a non-interacting Fermi gas. The notation for the 
AFQMC results is identical to Fig. [S] The diagrammatic MC 
results of Refs. [2ll . [2^ (solid up and down triangles) and the 
Bold Diagrammatic MC results of Ref. [2^ are shown as well 
(solid squares). The inset shows the vicinity of the superfluid 
phase transition at T^/ep ~ 0.15. 



to be a small effect for the Tan contact in Ref. [lO^, 
well as for the energy in Ref. 
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In spite of these shortcomings, the introduction of 
HMC into the AFQMC study of the UFG has largely 
solved the issue of sufficiently large spatial lattice dimen- 
sion and particle number N ^ which in turn has allowed 
calculations with a large particle number at lower densi- 
ties. DHMC studies for > 16 are in progress. These 
improvements will also apply to calculations away from 
the unitary limit. Finally, we would like to stress that the 
AFQMC method is entirely ab initio: once the coupling 
g is fixed by solving the two-body problem, no tuning 
with respect to experiment is required. While a more 
sophisticated analysis of the systematic errors cannot be 
provided at this point in time, the fact that theory and 
experiment agree reasonably well for both E/E^q and 
/i/e^ over a wide range of temperatures is both remark- 
able and encouraging. 
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